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Abstract 

We prove that, for m > 7, scalar evolution equations of the form ut = F(x, t,u, . . . , Um) which admit 
a nontrivial conserved density of order m + 1 are linear in Um- The existence of such conserved densities is 
a necesary condition for integrability in the sense of admitting a formal symmetry, hence integrable scalar 
evolution equations of order m > 7 are quasilinear. 



1. Introduction 

It has been shown that scale invariant polynomial of evolution equations of order greater than seven are 
symmetries of third and fifth order equations [1] and similar results are obtained in the case where negative 
powers are involved [2]. The problem of classification of arbitrary evolution equations is thus reduced to 
proving that such equations have desired polinomiality and scaling properties. As a first step in this direction, 
we study scalar evolution equation of general type and show that integrable equations of the form 

Ut^ F{x,t,U,. . . ,Um) (1-1) 

where F has partial derivatives of all orders, is linear in Um for m > 7. The integrability criteria that we 
are using is the one given in [3], namely the existence of a "formal symmetry" R, which is a formal series in 
inverse povers of D ^ d/dx, satisfying the operator equation 

Rt + [R,F4^0, (1.2) 

where F^ is the Frechet derivative of F. The solvability of the coefficients of R in the class of local functions 
requires that certain quantities denoted as p*-'-* be conserved densities. These conserved density conditions 
give overdetermined systems of partial differential equations for F and lead to a classification. 
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For m = 3, it is known that there are essentially nonlinear equations that are candidates for being 
integrable [3] and their study require special techniques as explored in [4] . Our proof based on the existence 
of the conserved density p^^^ is not applicable to third order equations, because for m = 3, the expression of 

given by (2.11) is no longer valid. For m > 7 the requirement that pi be a conserved density leads to 
two equations involving dF/dum which are are compatible only when F is linear im Um- For rn = 5 there 
are three such equations given by (4.16a-c), but these arc compatible and we cannot exclude the possibility 
of the existence of essentially nonlinear equations at the fifth order. 

In Section 2, we give a concise overview of the formal symmetry method and present the expression of 
the conserved densities p^'^ i = 1,2,3. In the proof of out quasilinearity result, we actually do not need 
the explicit expression of these conserved densities, because the only information we use is the existence of 
a conserved density of order m + 1 which is quadratic in Um+i- The computation of the top order terms in 
derivatives are given in Section 3. Section 4 is devoted to general results and to the discussion of third and 
fifth order equations. The main result is presented in Section 5. 

2. Notation and terninology, formal symmetries, conserved densities 

Let u = u{x,t). A function (p of x, t, u and the derivatives of u up to a fixed but finite order will be 
called a "differential function" [5] and denoted by ip[u]. We shall assume that ip has partial derivatives of all 
orders. For notational convenience, we shall denote indices by subscripts or superscripts in paranthesis such 
as in aj-j) or p^*) and reserve subscripts without parantheses for partial derivatives, i.e., 

, , d'^u du du 

u=u{x,t), uk = -^, = ^ ' 

, , dip dixi dw , 

ip=ip{x,t,u,ui,...,Un), ^''^d^' ^''^'dx' ^*^~di' ^ ' 

In Section 5 only, we shall use the notation p" = Pm+i,m+i) in order to simplify the presentation. 

If is a differential function, the total derivative with respect to x is denoted by Dip and it is given by 

n 

Dip = '^(Pi Ui+i + ip^. (2.2) 

1=0 

Higher derivatives can be computed by appling the binomial formula as 

n k—1 

= ('7')^'^^ + ^'"V.- (2.3) 

j=0 j=0 



Note that even if ip has an arbitrary functional form, D'^ip is polynomial in Un+i for i > 1. We shall denote 
generic functions that depend on at most m„ by 0{n). More precisely, 

ip = 0{un) if and only if dip/dun+k = for fc > 1 (2.4) 
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Clearly, if = 0(u„), then D^p is linear in and D'^if is polynomial in Un+i for i > 1. In certain places 
we shall need to know not only the order of a function but whether it is a polynomial or not. When this 
distinction is important, we shall use the notation y> = P{un), i.e., 

ip = P{un) if and only if ip = 0{un) and d^ip/du^ = Q for some fc. (2.5) 

For ^p and u as above, the total derivative with respect to t is denoted by Dt^p and is given by 

n 

Dt^^J2^iD^F + ipt. (2.6) 
Equalities up to total derivatives with respect to x will be denoted by =, i.e., 

ip^ip if and only if ip = + Drj, (2.7) 

for some differential fimction 77. 

A symmetry a of Ui = F[v] is a differential function which satisfies the linearized equation at = -F*"' 
where F^ is the Frcchct derivative of F while a conserved density p is a differential function satisfying pt = Drj 
for some differential function rj. The order of a symmetry or of a conserved density is the order of the highest 
derivative of u in cr or p. The recursion operator R is defined to be a linear operator that sends symmetries 
to symmetries, i.e. Ra is a symmetry whenever cr is a symmetry [5] and satisfies Eq.(1.2). A formal series 
which satisfies Eq.(1.2) up to a finite order is called a formal symmetry. 

In 1 + 1 dimensions the integral terms in the recursion operator can be expanded in inverse powers of 
D and one can look for a formal series satisfying Eq.(1.2). It is known that whenever i? is a solution, i?*^/" 
is also a solution of Eq.(1.2), hence the order of R is irrelevant as long as we deal with recursion operators 
expressed clS cL formal scries in Z? . 

It is known that the quantities 

^(-1) ^ p-l/m^ ^(0) ^ Frn_^IF^^ (2.8) 



where 



dF ^ dF 

Fm-i = (2.9) 



dUm Um-1 

are conserved densities for equations of any order [3] . 

For the computation of the conserved densities for arbitrary m, one can start with a formal series for R, 
substitute in Eq.(1.2) and solve the operator equation recursively, obtaining at each step linear first order 
differential equations for the coefficients of R. The solvability of these equations in the class of local functions 
gives the conserved density conditions, p[^^ = Dr]^''^ , including the ones above. 
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Alternatively, one can use the fact that the coefficient of in the commutator of any two formal 
series is a total derivative [3]. If we start with a formal series R of order 1 satisfying (1.2), the coefficient of 
will be the conserved density p^^\ By the remark above, R'' also satisfies (1.2), hence the coefficient of 
in R'' will also be a conserved density that we denote by p^^'h For k = m the top m — 2 terms in R^ 
should coincide with , and the coefficients of R are expressed in terms of the partial derivatives of F for 
each m. 

For third order equations, the expression of p^^^ involves rj'^^^^ defined by Dtp'^^^^ = Di]^-^^\ But for 
equations of order 5 or higher, p^^^ and p^^^ are independent of 77^"^^ and 77^°^. Similarly for 5th order 
equations, p^^^ depends on ri'^^\ but for m > 7, p^^^ is independent of r/^'^ with i < 2. We have computed 
the conserved densities p^^^ and p^^^ for m > 5 and p^^^ for m > 7 with the method described above and 
compared with the results of direct integration of the first order equations obtained recursively from Eq.(1.2). 
The formulas below were obtained analytically for m > 13 and proved by induction. We have computed 
the conserved densities at lower orders directly using REDUCE and a symbolic integration package [6] and 
checked that they obey the general formula. The results are given in Eqs.(2.11-13) where we use the notation 



(i) 



1/ -^rr 



1,2,3,4. 



(2.10) 



Then for m > 5 we have 



p(i) =a-\Daf 



i{m+\) 



Da a(i) + a 



24 



p(2) =a{Da) iDa(i) + ^a^^^ - j^^a^2) 



(2.11) 



(2.12) 



while for m > 7, 



+ 30a {Daf ^^^^^;§^Da^i) + ^^(i+i) ^^) - m(r^-i -1) ^(2) 



m(m^ — l){m+3) 



Da 



(m— l)(m— 3) 



a(i) + (m - 3)Da(2) 



(m-l)(2m-3) 3 , 6(m-2) r. 



+ 



60 . 
m(TO^-l)(m+3)' 



{Dan)f - ^Dan) a(2) + 



(3) 



(m-l)(2m-3) 



«(1) 



. (2m— 3) 2 I 8 I 4 2 8 

■ + m"(i)"(3) + ^"(2) - ^;;^a(4) 



(2.13) 
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3. Top four terms in the derivatives of differential functions 

In this section we shall obtain the expression of D'^ip up to top four highest derivatives. Let <^ and u be 
as in Eq.(2.1), and D''^p be given by (2.3). We start by writing the expressions of D''^p for A; < 6 up to top 
four derivatives. As there is a truncation in the summations, for fc > 4 we should assume that n > 3. 



Dip =^^ipi Ui+i + ipx, (3.1a) 

n n 

D'^ip =^ipi Ui+2 + ^Dipi Ui+i+ Dipa,, (3.16) 

i=0 i=0 

n n n 

D^(p = '^(pi Ui+3 + 2 ^ Dipi Ui+2 + ^ D'^ipi Ui+i + D'^tpx, (3.1c) 

j=0 i=0 i=0 

n n n n 

D^(p= ^ (pi Ui+4 + 3y^Dipi Ui+3 + 'Sy^D^ipi Ui+2 + y^^D^Vi ut+i + D^ifix + 0{un), (3. Id) 

i=n-3 i=0 i=0 i=0 

n n n n 

i=n—S i=n—2 i=0 i=0 

n 

+ ^Vi Wi+i + D^^x + OK+i), (3.1e) 

n n n n 

i=n—3 i—n — 2 i=n—l i=0 

n n 

^ Ui+i + £>Vx + 0(«n+2). (3.1/) 

i=0 i=0 

For k > 7, we keep terms of orders Wn+fe, Wn+fe-i, Mn+fe-2 and Un+k-3 and write the general formula as 
below. 

D if =ipn Un+k + fn-l Wn+fc-l + <Pn-2 Wn+fc-2 + fn-S Un+k-3 + ■ ■ ■ 
+ (fc - 1) [Dipn Un+k-1 + D(pn-1 Un+k-2 + D(p 
+ (''2 ^) [^^'Pn Un+k-2 + D^^n-l Un+k-3 + • • •] 
+ ^) [D^^pn Un+k-3 + ■■] 
+ 

+ (''2^) [D''~^^n Un+3 + D'^-'^tpn-l W„+2 + -D''"V«-2 Un+1 + • . ■ + D''-'^ ipi U4 + O'^-^ipo U3] 
+ (fc - 1) [D''-^iPn Un+2 + D''-^iPn-l Un+1 + D^-'^ipn-2 «„ + ... + D^-'^ipi U3 + D^-'^ipo U2] 
+ Vn Un+1 + Vn-1 + Vn-2 + . . . + D'^" Vl "2 + -D''" Vo ^i] 

+ D^-^^^+0{Un+k-A). (3.2) 

Note that, for k>7, only the indicated terms of the first four lines contribute to the top four terms. In the 
last three lines, there are contributions from top three terms in D''~^(pi, top two terms in D''~'^ipi and from 



the top term in D''~^ipi. For k>7 these different types of contributions do not mix up and we can obtain a 
general formula by inspection, obtaining recursively the top first, second and third terms of any derivative. 
Once the general formula is "guessed" it can be proved easily by induction. 

Proposition 3.1. Let (p and u be as in Eq. (2.1 ) and assume that n>3. Then 

D^ip =LPn Un+k + P{un+k-i), k > 2, (3.3a) 

D''(p =(pn Un+k + bn-1 + kD(p„]Un+k-l + P{Un+k-2), k > 3, (3.36) 
D''(p =ipn Un+k + [fn-l + kDip 

+ yn-2 + k Dipn-1 + {^^D'^LPn]Un+k-2 + P{Un+k-d,) k > 5, (3.3c) 
nj'^n+fc— 1 

+ [<Pn-2 + k D(pn-l + {'^)D'^ipn]u„+k-2 

+ [iPn-3 + k Dipn-2 + {2)D^<Pn-l + {3)D^V'n]Un+k-3 + P{Un+k-4), k>7 (3.3d) 

and tie expressions of D'^tp for k = 1, ... 6 are 

Dip =ipn Un+l + 0{un), (3.4a) 
D'^P =p>n Un+2 + [V'n-l + 2£)(p„]u„+i - U^+j + 0(w„), (3.46) 

D^(p =(pn Un+3 + + 3D</?„]m„+2 + [(pn-2 + 'dDipn-i + 'dD^(pn]u„+i 

- 3<^„„ Un+l Un+2 " 3£><^„„ ^n+l + Vnnn 

- 3v„,„_i u^+i + 0{un), (3.4c) 

D'^ip =(pn Un+4, + bn-1 + ADipn]Un+3 + ['Pn-2 + 4£>9?„_i + 6£)^(/J„]u„+2 

+ [(pn-3 + 4:D(pn-2 + QD'^'Pn-1 + 4-D^¥'„]u„+l 

- 4:ipnn Un+l Un+3 - ^D^Vnn - 12-Dv„„ "n+l Un+2 - 3ipnn "^+2 
+ ^Vnnn Un+2 - 10y„,„_i Un+l Un+2 - 12Dipn,n-l 

-t- ADipnnn ^n+l ^nnnn ^ ^pnn.n—l ^ri-j-l 

- ^'Pn,n-2 Un+l " ^^n-l,n-l M^+l + 0(u„), (3.4d) 
D^(p =(pn Un+5 + [i^n-1 + ^DiPn]Un+i + ['Pn-2 + 5£>y„_i + 10£>^y„]u„+3 

+ [pn-3 + ^Dipn-2 + 10£>Vn-l + 

nn Un+2 Un+3 ^ ^^DPnn Un+2 ~ 10'/'n,n — 1 ~^ P{Un+l)j (3.4e) 

D^p =p>n Un+6 + [y^n-1 + QDfn\Un+5 + ['Pn-2 + GDpin-l + 15£>^¥'„]u„+4 
+ [pn-3 + QDipn-2 + 15Z)^(^„_i + 2QD^ipn]Un+3 
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(3.4/) 



Proof. The expressions (3.4a-f) are obtained iteratively by expanding the derivatives in (3.1a-f) and keeping 
the relevant terms, which is a tedious but straightforward computation. 

Similarly, to prove (3.3d) by induction, one has to check that it holds for fc = 7, by replacing fc = 7 in 
(3.2). Then assuming (3.3d) holds for k >7, we compute 

ij'^+V ='fin Un+k+1 + Dip 

+ [Vn-l + k D(Pn]Un+k + [Difn-l + k D'^ ipn]Un+k-l 
+ [fn-2 + k D(pn-1 + {'^)D'^ipn]Un+k-l 

+ [Dipn-2 + k D'^ipn-l + {^^D'^ifn]Un+k-2 

+ [pn-Z + k Dipn_2 + (2)£'Vn-l + (^)£'Vn]Wn+fc-2 + P(«„+fc-3), (3.5o) 
=iPn Un+k + + (fc + 1) Dlf 



+ 
+ 



iPn-2 + {k + 1) D<Pn-l + (2) + ^ D"^ I Un+k-\ 
ip„-3 + {k + l)k D(pn-2 + Q+k D^(pn- 
+ P(u„+fe_3). 



(3) + (2) 



D^ifin 



Un+k-2 



(3.56) 



Using the relations 



k\ , (k + 1 
2 +^= 2 



k + l 
3 



in (3.5b), we obtain Eq.(3.3d). The validity of the remaining expressions in (3.3a-d) for the indicated values 
of k can be checked from Eqs.(3.4a-f). • 

4. General Results on Classification 

In this section we shall obtain certain general results on classification. We start by stating the following 
well known result which is proved easily from the conserved density condition. 

Proposition 4.1. Let F = F{x, t,u, . . . , U2k)- If p = p{x, t,u,..., m„) is a conserved density of the evolution 
equation ut = F, then pnn = 0, for n > 2 and A; > 1. 

Proof. Let n = 2k + I. Then Dtp is the sum of terms p2k+i-jD^''~^^~'' F ■ Integrating this k — j times 
we obtain D''~^ p2k+i-jD''~^^F which is quadratic in u^k+l-j- Thus the contribution to the highest order 
nonlinearity comes from the first term corresponding to Z = only. It follows that p„„ = 0. • 

For m odd, top 2 terms in the expansion of the derivatives contribute to the highest order nonlinearity 
and we have a nontrivial result. 
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Proposition 4.2. Let p be a conserved density of order n> m for an evolution of order m = 2fc + 1. Then 
for k>2, 

7TL 1 

— Fm Dpnn - {u - -m) DFm Pnn = PnnFm-l m>7, U > 171 (4.1) 

Proof. From the proof of Proposition 5.1 it will be seen that only top terms in Dip will contribute to the 
highest order nonlinearity. Let n = 2k + I + 1. Then the highest nonlinearity will be of order 3k + I + 1 and 

Dtp ^p„£>2fc+;+i^ ^ Pn-iD^'+'F + 0{ulk^i) 

[D'^+Vn - D'^pn-i] D^+'F + 0(mL+/). 

Furthermore, only top two derivatives of each term will contribute to the highest order nonlinearity. Hence 
we can use the formula (3.3b) which is valid beyond the third derivative. We note that (3.3b) can be used 
for fc = 2 also, because only the top derivative in D'^pn-i is needed. Thus 

(-1)*"+"^ Dtp = PnnU3k+l+2 + [k + l)DpnnU3k+l + l 

X FmU3k+l+l + [Fm-1 + {k + l)DFm\uzk+l "sfe+f+l + ^^C^Sfe+f) 

- iZ)(p„„F„) - Pnn[Fm-i + {k + l)DFj[ + (fc + l)F„£)p„„ wifc+i+i + 0{ulk+i) 

- [(fc + h)DpnnFm - {k + I + \)p„nDFm - p„nFm-l\ Uafe+J+l + <^('"3fe+;)- 

The coefficient of u\f._^i_^-^ gives (4.1) and the proposition is proved. •. 

The implications of (4.1) are listed below. Conserved densities of order n > m are necessarily quadratic in 
the highest derivative and there is essentially only one conserved density at each order. The corollary can 
be proved easily by replacing the indicated values for n in (4.1). 

Corollary 4.3. Let p be a conserved density of order n>m for an evolution equation of order m = 2k + l > 
5. Then 
i. For n = m, 

Fm Pmmm ~ -^mm Pmrri' (4.2ct) 

a. If p and rj are two conserved densities of order n with pnn 0, rjnn 7^ 0, then 

Pnn Vnn 

Hi. For n> m, 

Pnnn = 0. (4-2c) 

Applying (4.2a) to p^~^^ , we will show that there are three classes for F, the first one being the quasilinear 
equations. For m = 3, both the essentially nonlinear classes are well known candidates for integrable 
equations. 



Proposition 4.4. Let p = Fm^^™ be a conserved density of the evolution equation ut = F. Then F belongs 
to one of the classes below. 

i. pm=0- F = AUm + B, Am = Bm = 0, (4.3a) 

n. Pmm = 0- F={Aum + B)^-"' + C, Am=Bm = Cm = 0, (4.36) 

in- Pmm^O: Fm= [c(1)f2 + c(2)F + C^^)] "/^""'^ ^ ^(l) = ^(2) ^ ^(3) ^ q ^4 3^-, 

Proof. First note that Pm = ^ corresponds to Fmm = 0, in which case the conserved density condition 
(4.2a) is identically satisfied and gives the class (4.3a). Then we compute 

Pmm = ^(^ + 1)-Fm'" ^mm ~ m'^"^"' ^mmm (4.4) 

and pmm = gives 

(i + l)Fl^ + FmFmmm = (4.5) 

which can be integrated to give the class of equations (4.3b). For pmm 7^ 0, Eq.(4.2a) gives 

Fynm Pmmm 



Fm Pmm 



(4.6) 



which implies that F^ is proportional to Pmm and using Eq.(4.4) we obtain a third order ordinary differential 
equation. Integrating this twice we obtain 

Pm = P'F + V, IJ.m = t^m = 0. (4.7) 

Substituting p = F^^"^ in (4.7), integrating once more and renaming integration constants we have 
Eq.(4.3c). We also note that the class (4.3c) with c^^^ = corresponds to the class (4.3b). • 

The equation (4.3c) can be integrated for m = 3 but it has no antiderivative for m > 5. Thus for m = 3 we 
have the following well known result. 

Corollary 4.5 (Third order equations). Let Ut = F{x, t,u,..., U3). If p = F^^^^ is a conserved density 
then F belongs to one of the classes below. 

F = Aus + B, As = B3= 0, (4.8a) 

F = {Au3 + B)-^ + C, ^3 = B3 - C3 = 0, (4.86) 

^= iAui'+Bul+Cyf^ ^''^ ^3=i?3 = C3=i.3 = 0. (4.8c) 
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We now arrive at the discussion of equations of order m > 5. Recall that the expression of p^^^ given 
by Eq.(2.11) which is valid for m > 5 is of the form (4.2c) with 

Pnn = 2a- ^a^, (4.9) 

where n = m+l. Prom (4.1) we obtain 

flm 2 a 

which gives 

a = {Aum + B)-^ (4.11) 

where A and B are independent of Um- It follows that 

^^11 + m^^^^ ^- = i?m = C„ = 0. (4.12) 
Al — m 

It can be checked that this corresponds to 

ci = (t^^)'^', C2 = -2ciC, C3 = ciC^ (4.13) 

in (4.3c) and it is inconsistent with (4.3b). As we shall see that for m > 7 only quasilinear equations are 
allowed, we state the Proposition below for m = 5. 

Proposition 4.6 (Fifth order equations). Let p = p'-^^ given by Eq.(2.11 ) be a conserved density of the 
evolution equation Ut = F{x, t,u, . . . , u^). Then integrable equations consist of two classes 

F = Au5 + B, As = Bs = 0, (4.14a) 

i^ = -^(Awi+B)-2/3 + c, yl5 = B5-C5=0, A4 = Ci = 0. (4.146) 

Proof. Substituting (4.12) in the conserved density condition and taking logarithmic derivatives we obtain 

Fm-i/Frr, = mDA/A. (4.15) 

Taking the derivative of F with respect to Um-i and substituting, we see that the expression above should 
be linear in Um- The second derivative with respect to Um gives Cm-i = 0, while the coefficient of the linear 
term leads to Am-i = 0. Substituting m = 5, we have (4.14b). • 

Remark 4.7. Let p be a conserved density of order 6 for an evolution equation of order 5. In the expression 
Dtp, the coefficients of UgUQ, u^uq and u^Uq depend only on pee- These expressions give respectively 

7-^55 P66 — 5p665 = 0, 

— 42i^555 P66 — 2F55 P665 + ^^F^ P6655 = 0, 

— 28^5555 P66 + 24i^555 P665 — 23i^55 P6655 + 15^5 P66555 = 0. 
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After eliminations we obtain 

-F55 P66 =1^5 P665 = 0, (4.16a) 

^55 P665 =5-F5 P6655 = 0, (4.166) 

-P55 P6655 =^-^5 P66555 = 0. (4.16c) 

It can be checked that the equations above are consistent with 

where the last equation corresponds to Eq.(4.10) with m = 5. Thus we cannot exclude the possibility of the 
existence of essentially nonlinear equations at the 5th order 

5. Equations of order greater th£in seven. 

We consider now evolution equations of order m > 7 admitting a conserved density p of order n = m + 1, 

i.e., 

m = 2fc+l, n = m + l = 2k + 2. 
From Corollary 4.3, we know that p is quadratic in «„, i.e, 

P = aul + pUn +7, Q!n = = 7n = 0. (5.1) 

For convenience, in the proof of Lemma 5.3 and afterwards, we shall use the notation 

p' = Pm+l, P" = Pm+l,m.+l, Pj=Pm+l,j, Pj = Pm+l,m+l,j, j = 0,...m. 

The following relation can be easily deduced from Eq.(5.1). 

Lemma 5.1. Let p = p{x, t,u,..., Um+i) be a quadratic function of Um+i- Then 

Dp' =p'^ Um+l + 0{Um), D'^P" = Pm Um+2 + P(Um+l), p" = Pm + P{Um+2), 

Dp'j =p'^ Um+2 + P{Um+l), D'^ p'j = p'^ Um+3 + f (Um+2). (5.2) 

We shall now show that, in the integral of Dtp, the contribution to top two nonlinear terms u'^k_^2 ^'^^ 
"sfc+i ^^^^ come from top four derivatives. 

Proposition 5.2. Let p = p{x, t,u,... u„) and Ut = F{x, t,u,..., Um), where n = m + 1 and m = 2fc + 1. 
Then 

(-1)'=+! Ap = [d^+^Pu - -DVn-i] D^^^F - D^F - D'^- Vn-s] D'^ F + 0{usk)- (5.3) 
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Proof. We write 



Dtp =pn D^'F + pn-i D^'-^F + pn-2 D""-^ F + p„_3 D""-^ F + pn-4 D^-^F + ...+poF + pt, 
=P2k+2 D^^+^F + p2fe+i D-'^+^F + p2k D^'^F + p2k-i D^'^-^F + p2fe-2 D^'^-^F + . . . 
+ PqF + pt, 

We integrate the first term k+1 times, the second and the third k times etc., so that each term is a product 
of terms of orders equal or differing at most by one. Hence we have 



{-l)''+'Dtp ^ D'^+'p2k+2 - D''p2k+1 



D'^+'^F - 



D''p2k - D''-^p2k-i 



D^'F 



- D^-^p2k-2-D''-'p2k-s D^-'F + .-. + poF + pt. 



lfe-2 



It will be sufficient to show that third term is of order u^ki because by a similar reasoning it can be seen 
that the remaining terms are of lower orders. As p2k-2 is of polynomial in U2k+2 the term in the brackets is 
linear in u^k+i and polynomial in other derivatives of order larger than m. Similarly D''~^F is polynomial 
in usk, hence after integration by parts, their product is of order usk- It follows that the contribution to top 
two nonlinear terms comes from top four terms in the expansion of Dtp. • 

The contribution to the coefficient of u'^f^_^_2 comes from the top two derivatives of the first term of (5.3) 
and it has already been given in (4.1). We shall need now the coefficients of 

■"3fe+2 U2k+2 and U3fe+1 U2k+4- 

We shall now show that at most top four terms of each derivative contribute to the coefficient of u'^k+i ^^'^ 
actually the contribution to the coefficient of u1i._^-^^U2k+4 comes only from the top two terms in each product 
in 

{D'^+'pn - D^Pn-l)D^+^F 

For this we shall compute the relevant terms in each factor of Dtp separetely. 

Lemma 5.3. Let pOp{x,t,u, . . . ,Un) be quadratic in Un = Um+i, F = F{x,t,u, . . . ,Um) and let k > 6. 
Then 

D^+'^p' - D^Pm =p" Um+k+2 + [r\ p'L Wm+l + 0(Wm)] Um+fe+1 

+ [r2 Mm+2 + P(Mm+l)] '^m+k 

+ [^3 Mm+3 + ^'(Wm+2)] Wm+fe-1 

+ P{Um+k-2) (5.4a) 
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where 



D'^Pm-l - Vm-2 =P{Um+l) Wm+fe+1 + P{Um+2) Um+k + P{Um+k-l) (5.46) 

Um+\ + 0{u 
+ [/2 Um+2 + P{Um+l)] 

+ [/s Fmm Um+3 + P{Um+2)] '^m+k—2 

+ P{Um+k-3) (5.4c) 

D'^F =0{Um) + P{Um+l) Um+k-1 + P{Um+k-2) (5.4d) 

r,=k + l, r, = l+(^f), rs = (^r) + ( - (2) (5-5a) 
/i=fc + l, /2=('r). /3 = (T) (5-5^>) 
Proof. We substitute if = p' in (3.3d) which is vahd for A; + 1 > 7 and get 

V =P" Um+k+2 + [p'm + {k + l)Dp"]Um+k+l 

+ [p'm-1 + ik + l)Dp'm + ('t')^V"] Wfc 

+ [p'm-2 + l)Dp'^_, + {'t')D^p'm + ('t')l^V"] Um+k-1 

+ P{Um+k-2) 

=p" Um+k+2 + [p'm + + l)PmUm+l + 0{Um)]Um+k+l 



+ 



[(fc + i) + ('=+i)]p>™+2 + p(u™+i; 

[{'X') + {'f)\p'Lura+Z + P{Um+2) 



Um+k 
Um+k-1 



+ P{Um+k-2). (5.6) 

Then we repeat the same computations for D'^pm up to top three terms using (3.3c) which is valid for A; > 5. 

D'^Pm =Pm Um+k+1 + [Pmm + kDp'^]Um+k 

+ [Pm,m-1 + kDpmm + {2)D'^P'm] Um+k-1 
+ P{Um+k-2) 
=P'm Um+k+1 + [kpmUm + 2 + PiUm+l)]Um+k 



+ 



Um+k-1 



{2) PmUm+3 + P{Um+2) 
+ P{Um+k-2). (5.7) 

Subtructing (5.7) from (5.6) we obtain (5.4a). The derivation of Eqs.(5.4b-d) is straightforward. • 

In the products of (5.4a) with (5.4c) and (5.4b) with (5.4d) most of the terms will not contribute to the 
top two nonlinear terms, as indicated below. 
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Lemma 5.4. For k > 6, the following relations hold. 

Um+k+2 P{Um+k-3) = U^+k+l P{Um+k-2) = Um+k P{Um+k-l) = P{Um+k-l) (5.8a) 
Um+k+2 Um+k-l P{Um+l) = Um+k+2 Um+k-2 P{Um+2) 

— Um+k+l Um+k P{Um+l) = Mm+fc+1 Um+k-l P{Um+2) = "m+fe P{Um+2) (5.86) 

Um+k+2 Um+k V (Wm) — Mm+fe+l 0{Um) (5.8c) 

It follows that the contribution to the top two nonlinear terms of Dtp comes from 

P" Um+k+2 + Pm Vl Um+1 Um+k+l + r2 Um+2 Um+k + r3 Um+3 Um+k-l] 

Fm Um+k+\ + Fmm [/l Um+1 Um+k + /2 Um+2 Um+k-l + /s Um+3 Um+k-2] , (5.9) 

where r,'s and /j's are given by Eqs.(5.5a-b). Substituting these values and integrating by parts we obtain 
the following result. 

Proposition 5.5. Let Ut = F(x, t,u, . . . , Um) be an evolution equation of order m and p = p(x, t,u, . . . , u„) 
be a conserved density of order n = m+1 satisfying Pnnn = 0. The coefRcients of M3^_|_2M2fe+2 and u^f._^_-^U2k+'i 
give respectively 

{2k + l)p':^Fm = {2k + i)p"Fmm, (5.10a) 
{2k + l)(fc2 + k + 6)p'^Fm = {2k + 3)(fc + l)(fc + 2)p"Fmm (5.106) 

where p" = pnn 

Proof. The general formulas can be used for fc > 6, i.e., m > 13 and Eqs.(5.10a-b) can be obtained 
from Eq.(5.9). For m = 5, straightforward computation of the conserved density condition is possible with 

REDUCE and the results have already been given by Eqs.(4.161-c). For m = 7 and m = it is necessary to 
use a combination of symbolic and analytic computations. For m = 11. there is little discrepency from the 
general formulas and the required coefficient scan be obtained easily with analytical. • 

It can be easily checked that the equations (5.10a) and (5.10b) are inconsistent except for k = 2. For 
k = 2, (5.10b) is the coefficient of u^us, and leads to Eq.(4.16b) after integration. Thus for fc > 3 Thus 

Q^m Fm — Oi Fmm — 0. 



It follows that for a ^ 0, Fmm = 0, and we have the following corollary. 
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Corollary 5.6. Let Ut = F{x, t,u,..., Um) be an evolution equation of order m>7 and p = p{x, t,u,..., u„) 
be a conserved density of order n = m + 1, with pnn 7^ 0. Then Fmm = 0. 

Note that in order to prove the quasihnearity resuh, the exphcit form of p'^-* is not needed. We have 
only used here the fact that it is indeed quadratic in Um+i- 
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